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The initial steps of MBE growth of GaAs on /32-reconstructed GaAs(OOl) are investigated by performing total energy 
and electronic structure calculations using density functional theory and a repeated slab model of the surface. We study 
the interaction and clustering of adsorbed Ga atoms and the adsorption of As2 molecules onto Ga atom clusters adsorbed 
on the surface. The stable nuclei consist of bound pairs of Ga adatoms, which originate either from dimerization or from 
an indirect interaction mediated through the substrate reconstruction. As2 adsorption is found to be strongly exothermic 
on sites with a square array of four Ga dangling bonds. Comparing two scenarios where the first As2 gets incorporated 
in the incomplete surface layer, or alternatively in a new added layer, we find the first scenario to be preferable. In 
summary, the calculations suggest that nucleation of a new atomic layer is most likely on top of those surface regions 
where a partial filling of trenches in the surface has occurred before. 



I. INTRODUCTION 

Growth of well-ordered crystals of III-V compound 
semiconductors requires the incorporation of both 
constituents in the correct stoichiometric amounts. 
For instance, in MBE growth of GaAs, the Ga atoms 
and As molecules provided by the beam sources must 
be adsorbed and incorporated into the growing sur- 
face in such a way that the surface stoichiometry is 
locally maintained. The atomistic processes behind 
stoichimetric growth are complex and not yet fully 
understood at present. Moreover, the substrate sur- 
faces used for growth of arsenide compound semicon- 
ductors, in particular the frequently used GaAs(OOl) 
surface, show a variety of complex surface reconstruc- 
tions. Under moderately arsenic-rich conditions, as 
are commonly used during growth, the GaAs(OOl) 
surface displays reconstructions with a (2 x 4) sym- 
metry: the a, (3 and p2 reconstructions, which con- 
tain surface As dimers as common building blocks. 
The strongly corrugated (32 reconstruction, which ex- 
poses three layers of atoms, prevails in a wide range 
of growth conditions and serves as the starting config- 
uration for growth on the GaAs(OOl) substrate. For a 
well-controlled growth, it is required that this struc- 
ture recovers after film deposition, at least after a 
short growth interruption. However, it was already 
understood in early growth models [jlj that different 
reconstructions may appear locally on the growing 
surface, acting as metastable intermediates before a 
newly grown layer is completed. The details of these 
structural transformations remained unclear until re- 
cently. Only with the help of detailed STM studies 
Q has it become possible to refine our understand- 
ing of the elementary steps of growth ||. STM pic- 
tures taken from samples after submonolayer deposi- 
tion followed by a fast quench to room temperature 
show two major processes contributing to growth on 
the /32-reconstructed surface: the filling of trenches 
and the formation of small islands that later become 
part of the top-layer As dimers in the new layer. In 
the present paper, we present results of ab initio cal- 



culations for the atomistic steps of these two growth 
scenarios. 

II. CALCULATIONS 

We performed ab initio calculations using density 
functional theory to describe exchange and correla- 
tion in the electronic many-particle system. A recent 
version [|| of the generalized-gradient approximation 
for the exchange-correlation functional was employed. 
All calculations were done with slabs consisting of 
seven or eight atomic layers and a 4 x 4 lateral unit 
cell. The bottom layer was passivated with pseudo- 
hydrogen atoms and kept fixed, while the top six or 
seven layers and adatoms were allowed to relax. 

Our calculations use fully separable, norm- 
conserving pseudopotentials |i|-[7j to describe the 
electron-ion interaction, constructed from an all- 
electron atomic calculation with the GGA functional 
[Q. Gonze's analysis § was used to confirm that un- 
physical ghost states were not present in the separa- 
ble representation. The wave functions were expanded 
in a plane wave basis || with a cutoff energy of 10 
Ry, and the k-space integration was performed with 
a special k-point set, with a density equivalent to 64 
k-points in the Brillouin zone of the (lxl) surface 
unit cell. 



III. ADSORPTION 

Both in thermodynamic equilibrium and over a 
wide range of growth conditions, the GaAs(OOl) sur- 
face forms the (32(2 x 4) reconstruction. Each (2 x 4) 
unit cell is built up from two As dimers and two miss- 
ing dimers in the topmost layer and a missing Ga pair 
in the second layer (see Fig. |l]). The missing atoms 
give rise to trenches running in the [110] direction 
separated by mountains of adjacent As dimers in the 
top layer. 
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FIG. 1. Adsorption sites for Ga atoms on the /32(2 x 4) 
reconstructed GaAs(OOl) surface. The As and Ga atoms 
of the substrate are displayed as white and black circles. 
The unit cell of the reconstruction is marked by the dashed 
rectangle. Adsorption sites are labeled within the 4x4 cell 
used in the calculations. 
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The As atoms in the third layer exposed in the 
trenches also form dimers. On the atomic level, one 
can imagine two principal ways in which growth can 
proceed on this surface. One possibility is the nucle- 
ation of new layers on top of the existing 'mountains'. 
Alternatively, the trenches could be filled up first, ei- 
ther partially or completely, before nucleation of new 
layers starts afterwards in these surface regions. 

The adsorption of single gallium atoms on the 
GaAs(OOl) surface has been studied previously by 
means of density functional theory calculations |Tc[ |. 
On the /32-reconstructed surface, a Ga atom preferen- 
tially adsorbs between two As dimers in line with the 
dimer axis, at adsorption sites Ai in the trench or in 
A2 in the top layer, see Fig. [jjand Fig. |2|a. In these cal- 
culations, the substrate atoms were allowed to relax 
after deposition of the Ga adatom, but the bonding 
topology of the substrate atoms was maintained. A 
different kind of adsorption sites arises when break- 
ing of substrate bonds is taken into consideration. The 
Ga atom may split an As dimer in the trench (A3 , see 
Fig. [|b) or in the top layer (A4, see Fig. ||c) and ad- 
sorb in a two-fold coordinated site JToJ . Ga atoms in 
these sites are much more strongly bound than in the 
Ai and A2 sites. Adsorption there already constitutes 
the first step to incorporation of the Ga atoms. More- 
over, the binding energy is higher in the site A3 in the 
trench than in the A4 mountain site. As an alterna- 
tive to the adsorption sites Ai and A3, the Ga atom 



bonding to the As atoms in the trench may tilt away 
from the As dimer axis and form an additional bond 
with a Ga atom at the side wall of the trench (Bi 
and B3, see Fig.|d). In this way the Ga atom reaches 
a three-fold coordination. Our calculations with the 
GGA functional show that the three-fold site B3 ob- 
tained by splitting the As dimer is less strongly bound 
than the corresponding two-fold site A3, while a Ga 
atom between As dimers experiences a stronger bind- 
ing in the Bi site than in the Ai site (see Table |). 




c) ~ d)" 

FIG. 2. Schematic illustration of some adsorption sites 
of Ga and accompaning substrate rearrangements on the 
GaAs(OOl) /32 surface: a) Ga adatom in site Ai between As 
dimers in the trench, b) two- fold coordinated Ga adatom 
in site A3 splitting an As trench dimer, c) two-fold coor- 
dinated Ga adatom in site A4 splitting an As dimer in 
the top layer, d) three-fold coordinated Ga adatom in site 
B3 in the trench with additional bond to one side wall 
(dashed). As atoms are shown as white circles, Ga atoms 
as black circles. 
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IV. NUCLEATION 

MBE growth of GaAs is usually performed by ap- 
plying a flux of arsenic molecules that exceeds the 
flux of Ga atoms. Under these conditions the growth 
rate is mostly controlled by the diffusion and incor- 
poration of the Ga atoms, while As is easily avail- 
able everywhere on the surface. Therefore we focus 
on the interaction of Ga adatoms on the surface that 
is mainly responsible for nucleation. 

Gallium adsorption at site A3 has the largest bind- 
ing energy. Therefore, in the limit of low Ga coverage 
and low mobility, most of the deposited Ga atoms will 
get incorporated at randomly distributed A3 sites. At 
growth temperatures, however, thermally activated 
jumps of Ga adatoms from A 3 into neighboring sites 
occur frequently, with a rate of ~ 10 5 s _1 Jll]]. In this 
way, other adsorption sites, though being higher in 
energy (see Tab. |), will become populated according 
to a thermal distribution. If Ga adatoms in neighbor- 
ing adsorption sites interact attractively, clusters of 
Ga adatoms will form which act as precursors of is- 
land growth. Without substantial attraction, entropic 
effects lead to a preference for isolated single adsorbed 
Ga atoms in the small coverage limit. 

We have studied the interaction between Ga 
adatoms for various configurations. The interaction 
energies are defined with respect to the adsorption of 
individual Ga atoms in isolated A3 sites: 

AE{A 1 A N ) = E(A U ...A N )- NE(A 3 ) (1) 
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TABLE I. Interaction energy of clusters of Ga adatoms 
in the trenches (left) and in the top layer (right) of the /32 
reconstruction. All numbers are given relative to isolated 
Ga adatoms in A3 sites, with negative numbers indicating 
attractive interaction, with accuracy 50 meV. 
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FIG. 3. 'Trench' adsorption geometries (cf. left column 
in Tab. for Ga atoms in a trench: a) two Ga atoms in 
sites A3A3 splitting two As dimers, b) in adjacent B3B1 
sites, thereby splitting only one As dimer. 



The left column in Table | shows the interaction 
energy for Ga atoms adsorbed in the trench. The Ga 
atoms alternatingly occupy adjacent sites of the type 
Ai and A3, or Bi and B3. In this way, Ga-As-Ga 
...chains in [110] direction are formed from the As 
dimers in the trenches. From the calculated energet- 
ics, we derive the following growth scenario: When a 
diffusing Ga adatom approaches another Ga adatom 
in an A3 site, they can stabilize each other by forming 
a B3B1 pair. This process results in a local a recon- 
struction in one particular (2 x 4) unit cell. Although 
the Ga atoms do not form bonds with each other, 
such a pair is bound by 0.6 eV. Attachment of a third 
Ga adatom increases the interaction energy by only 
0.4 eV, to 1.0 eV. The unexpected stability of a Ga 
pair is due to the fact that an As dimer bond is re- 
stored when the diffusing Ga atom moves from the 
next-nearest neighbor site to the nearest-neighbor site 
with respect to a Ga atom adsorbed in the A3 site (see 
Fig. H for illustration) . 

Since the clustering of only two Ga atoms already 
results in a large increase in the binding energy, we 
expect the rate of formation of such pairs from a lat- 
tice gas of single diffusing Ga adatoms to exceed the 
rate of break-up of pairs even for a moderate super- 
saturation of Ga on the surface. Thus our calculations 
demonstrate that the Ga adatom pair in the trench 
acts as a stable nucleus in the sense of nucleation the- 
ory. 

Next we consider the alternative scenario, nucle- 
ation of a new layer without previous trench filling. 
The calculated results are collected in the right col- 
umn of Table |. In general, Ga adsorption in the top- 
layer As dimers is energetically less favorable than 
adsorption in the trench dimers. Therefore popula- 
tion of the top-layer sites is considerably smaller for 
low coverages, and only increases when the trench 
dimer sites are mostly occupied. However, when two 
Ga atoms adsorb in neighboring parallel As dimers, 
they reach a stable configuration due to the forma- 
tion of a Ga dimer (A4A 4 in Tab. Q). The decay of 
the Ga dimer into two single Ga adatoms in A3 sites 
is an endothermic process. This relative stability is 
maintained for those larger clusters that allow for for- 
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mation of Ga dimers oriented perpendicularly to the 
previous As dimers in the layer below, i.e. the A4A2A 4 
and A4A 4 A 4 A 4 clusters. We note that the latter clus- 
ter has twice the binding energy of a single Ga dimer, 
A 4 A 4 . Thus there is no extra attractive interaction 
between Ga dimers. If larger islands of Ga addimers 
form on surface regions where the trenches are already 
filled up with Ga atoms, this is not due to an attrac- 
tive interaction. However, such Ga islands may form 
accidentally and survive due to a limited mobility of 
the dimerized Ga atoms constituting these clusters. 

After a cluster of Ga adatoms has formed in the 
trench, any further Ga adatoms that are deposited on 
the adjacent mountains will be more likely to remain 
on the mountains long enough to form Ga dimers 
there, since they would have to move further to find 
empty trench sites to occupy. This suggests that the 
formation of a Ga cluster in the trench acts to pro- 
mote the subsequent formation of Ga dimer pairs on 
the adjacent mountains. 

In contrast to the situation in the trench, the 
formation of Ga-As-Ga . . . chains in [110] direction 
by splitting top-layer As dimers is associated with 
only a minor gain in binding energy (clusters A4A2, 
A 4 A2A 4 ). Relative to individual Ga atoms sitting in 
top- layer As dimers (A4 sites), forming a chain con- 
taining two (three) Ga atoms gives a binding energy 
of 0.2 eV (0.5 eV). However, these structures are un- 
stable against decay into individual Ga adatoms in A3 
trench sites, and therefore should have little relevance 
for growth. 

V. THE ROLE OF ARSENIC 

Up to now, we have not considered the possibility 
of enhanced stability of the above structures due to 
adsorption of arsenic. As sources of arsenic in MBE 
growth, both AS2 and AS4 molecular beams are in 
use. For the issue of enhanced stability due to ar- 
senic, it is sufficient to consider the simpler case of 
AS2 adsorption. While AS4 is believed to split into 
fragments upon adsorption, the AS2 molecule can be- 
come incorporated into the (32(2 x 4) structure with- 
out dissociation. Ab initio calculations [ ^2|JT^ ] show 
that the binding energy of chemisorbed AS2 depends 
very much on the local environment. On an ideal (32- 
reconstructed surface, an arsenic addimer is bound 
to the top-layer As dimers by 1.65 eV JlJ]. It can be 
shown that As ad-dimers bound to Ga atoms will stay 
permanently adsorbed even at standard growth tem- 
peratures, while As dimers that bind to As atoms or 
to only one Ga atom are more weakly bound and will 
either desorb or react with diffusing Ga atoms |l2]j . 
This is consistent with the experimental observation 
that arsenic incorporation at these temperatures only 
proceeds when the surface is simultaneously exposed 



to a Ga beam providing excess Ga adatoms on the 
surface fl4|| . 

We expect that enhanced stability due to AS2 ad- 
sorption is most relevant for those structures contain- 
ing L-shaped patches of three Ga atoms or rectan- 
gular patches of four Ga atoms in adjacent sites. In 
these local environments, an adsorbing AS2 molecule 
can attach to the surface by building up three or 
four As-Ga backbonds. For the nucleation scenario 
in the trenches, such a situation occurs already for a 
single Ga atom adsorbed in the trench, for instance 
in the A3 site. The AS2 forms one backbond to this 
Ga atom, while simulataneously backbonds with two 
other dangling orbitals of Ga atoms at the sidewall 
of the trench are established. The binding energy of 
an AS2 molecule relative to the gas phase is 1.9 eV 
at this site, too low to make it permanently adsorbed 
at frequently used growth temperatures. When two 
Ga atoms are on adjacent sites in the trench, like in 
the local a structure (cluster B3B1, see also Fig. ||b), 
the adsorbing AS2 can establish four backbonds and 
will transform this structure to the (3 reconstruction, 
that contains three parallel As dimers (see Fig. [|c). 
The binding energy of AS2 is 2.4 eV in this environ- 
ment |l^]. Further AS2 adsorption on the local (3 re- 
construction, which would lead to a complete filling 
of the trench, has been found to be energetically un- 
favorable. We find a binding energy of only 0.9 eV for 
an As dimer filling in the fourth and last As dimer site 
in the top layer of the (2 x 4) unit cell. Since this bind- 
ing energy is much lower than the binding energies for 
an As addimer in other sites which we have found to 
be ultimately unstable at standard growth tempera- 
tures |l2j , we conclude that the filling in of the fourth 
dimer site on the local (3 surface should not play a 
major role as an intermediate configuration during 
standard growth of the (32 surface. Moreover, we find 
it to be more favorable for an AS2 molecule adsorbing 
on a local (3 structure to attach itself onto the As top 
layer, with its axis oriented perpendicularly to the ex- 
isting As dimers. We expect that complete filling of 
the trench does not occur until nucleation of the new 
mountains of the next layer up changes the structure 
so that it is no longer locally the (3 structure. 

We note that it may not be easy to distinguish ex- 
perimentally between a three-dimer (3 reconstruction 
and a four-dimer completely filled trench: In STM im- 
ages, the filled dangling bond orbitals of the As dimers 
in the top layer extend out far enough that the bright 
stripes corresponding to the two As dimers of the 
mountain appear considerably wider than the dark 
stripes corresponding to the two missing As dimers 
of the trench in the (32 regions. It seems quite plau- 
sible that the three dimer structure would appear as 
if it had a completely filled trench in STM pictures, 
since the filled dangling bond orbitals of the three As 
dimers in the top layer could extend far enough out 
to mask the fourth empty dimer site. 
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For nucleation of a new layer, we consider the pos- 
sibility that the three-Ga atom cluster A4A2A 4 and 
the four-Ga atom cluster A4A2A 4 A 2 (see Fig. [H|b) 
could gain in stability by getting "capped" with As2 
(Fig. |5|c). Similar to adsorption in the trench, we find 
that an As2 molecule with only three backbonds is 
rather weakly bond, by 1.7 eV, and thus can only play 
the role of an intermediary species in growth. How- 
ever, for the four-Ga atom cluster in the top layer, we 
find that As2 binds even more strongly there than on 
a cluster of Ga adatoms in the trench. The binding 
energy for an As addimcr on this four-atom cluster 
on the mountain is 2.7 eV. This indicates that un- 
der the usual growth conditions, any such four-atom 
cluster which forms is likely to be rapidly "capped" 
by an As addimer, becoming a very stable nucleus 
for the mountains of the next layer up. Since the Ga 
adatoms which fall on any local area are likely to mi- 
grate to the trench sites before such a four-atom clus- 
ter on the mountain has a chance to form, until the 
nearby Ga sites in the trench are completely filled, we 
expect that growth will generally proceed by a par- 
tial filling of the trenches, and formation of local re- 
gions of the (3 structure, followed by nucleation of the 
mountains of the next layer up in regions adjacent 
to locally filled trenches. This suggested growth se- 
quence is shown in Fig. ^|. Since nucleation of the new 
mountain is expected to proceed before the fourth As 
dimer adsorbs, completely filling the trench, we see 
that the new mountain must nucleate above the cen- 
ter or above the sloping sides of the mountains of the 
original layer - not above the center of the original 
trench. This may explain why the new mountains are 
observed to grow above the old mountains in STM 
pictures ||, instead of above the trenches. 

To make quantitative statements about the impor- 
tance of arsenic adsorption for stabilizing structures 
during growth, we need to take into account the actual 
conditions in the growth chamber. In the next sec- 
tion we discuss how this can be achieved in a growth 
model. 



a) 




FIG. 4. Growth scenario for filling of the trenches on 
the clean (32 surface (a), via formation of Ga atom pairs 
B3B1 (b), followed by As2 adsorption that leads to a local 
/3 reconstruction (c). The pictures represent show part of 
the slab in side view, with relaxed atomic geometries as 
obtained from the calculations. 
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c) 



FIG. 5. Growth scenario with nucleation of a new 
layer, starting from the clean f32 surface (a), nucleation 
of the four-atom Ga cluster A4A2A4 A 2 (b), and As 2 ad- 
sorption on the Ga cluster (c). 
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FIG. 6. Growth scenario proposed as a summary 
of the calculated energetics (schematic): the local /3 
reconstruction (b) acts as a precursor for nucleation 
of the new layer (c or d). 



VI. IMPLICATIONS FOR A GROWTH MODEL 

In empirical growth models for GaAs, often only 
the Ga atoms are treated explicitely [|15|. This ap- 
proach is justified if the concentration of As atoms 
on the growing surface is completely determined by 
thermodynamic equilibrium. In MBE growth the flux 
of incoming As molecules is much higher (typically 
by a factor of five) than the number of As atoms per 
unit area and time that are eventually incorporated 
into the crystal (which is equal to Ga flux). There- 
fore the net fluxes of adsorbing and desorbing As2 
molecules almost balance, and the concentration of 
As atoms can be considered as being close to equilib- 
rium |l6f| . In particular, this may be true if the res- 
idence time of the arsenic molecules is enhanced by 
a physisorbed precursor state fl4|| . Under the equilib- 
rium assumption, the chemical potential of arsenic, 
fj,As(p,T), is the same both inside the crystal and 
in the gas phase, and is determined by the sample 
temperature T and arsenic partial pressure p in the 
growth chamber. To be specific, we consider growth 
from As2- The gas phase chemical potential can be 
written down explicitely, 



Ma? (P,T) 



~&coh 



TS vib (T) 

TS Iot (T) - T(5 tl . a „s(p, T) - 5/2 k). (2) 



For the cohesive energy of As2, we employ our cal- 
culated result of 2.1 eV/atom, which is in fair agree- 
ment with the experimental value of 1.96 eV/atom 
p7[ . The entropies of vibration, rotation and trans- 
lation are calculated using the standard textbook ex- 
pressions [0 , with the experimental data on spectral 
properties as input data. 



Due to the strong binding of Ga atoms in A3 sites, 
the evaporation of Ga adatoms is negligible compared 
to the incoming flux for temperatures below 1000 K. 
In contrast to ref. (T^| , we therefore argue that the lo- 
cal concentration of Ga adatoms cannot be described 
by an equilibrium assumption, but should be treated 
by a kinetic model. The incoming flux of Ga atoms 
drives the system out of equilibrium. Only deep in 
the GaAs bulk, a chemical potential for Ga atoms, 
/1q| As bulk , can be defined by the law of mass action 
for formation of GaAs from the elements, 

^GaAs bulk^ T) = MGaAs(T ) _ MAs(p , T ). (3) 

Substituting ha s (p,T) from eq. ||, this equation de- 
termines the free enthalpy per particle a Ga atom 
finally reaches after incorporation. For various ad- 
sorption configurations, labelled i, the free enthalpy 
per particle g^ a is still well-defined even under non- 
equilibrium conditions, but no chemical potential can 
be specified. Differences between goa for Ga atoms in 
the beam, at the surface and in the bulk act as driving 
force for incorporation and growth. 

In Fig. 0, we compare ^Ba As bulk to ^Ga in the bulk 
of elemental gallium. By virtue of eq. ||, ^Ga^ bulk be- 
comes a function of temperature and AS2 background 
pressure. For the condensed phases of both GaAs and 
elemental Ga, we have assumed that the pressure de- 
pendence of [i may be neglected. The temperature de- 
pendence of \x for the condensed phase is calculated 
from a Debye model for the lattice vibrations, 



(4) 



MGaAs(T) = -2Elf 3 + 2 J dT' c G aAs (T')^^ (5) 

c v {T) = 9k B (T/T D ) 3 £° /T dx {e ff 1)2 - (6) 

Here -E C oh is the calculated cohesive energy per atom, 
3.1 eV for GaAs and 2.6 eV for Ga, and cy(T) is the 
specific heat per atom of GaAs and of the bulk phase 
of Ga, respectively. The Debye temperature To was 
taken to be 344 K for GaAs p| and 240 K for Ga 

& 

Fig. [7] shows that a background pressure of AS2 in 
the range of 10 -4 Pa to 10~ 3 Pa is required to stabilize 
the GaAs crystal under typical growth temperatures, 
between 700 K and 800 K. For 10~ 3 Pa As 2 pres- 
sure, the GaAs surface will become unstable against 
formation of gallium droplets at temperatures above 
- 800 K. 

The free enthalpy per atom introduced above allows 
to compare the energetics of structures on the kinetic 
pathway of growth that contain different amounts 
of gallium and arsenic. From our ab initio calcula- 
tions, we obtain the relative energies Ae^-* of vari- 
ous adsorbate structures on GaAs(OOl) at p — and 
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T = with respect to a reservoir of elemental Ga. 
The transfer of a Ga atom from the reservoir into 
the A3 site on the surface is endothermic by 0.55 
eV. The values of AeW can be read from Tab. | us- 
ing the relation Ae (i) = 0.55eV + AE (i ^ /N Ga , where 
iVca is the number of deposited Ga atoms. Formally, 
each deposited As atom goes along with the trans- 
fer of a Ga atom from the Ga bulk to the GaAs 
bulk reservoir, leading to a free enthalpy change equal 
to ^l Ashulk (p,T) -/z Ga (T). Hence we can directly 
compute the quantity 



3 «+A S «T = MGa (T) + Ae« 



GaAs bulk 



(p,T)- MGa (T))7V As /iV, 



Ga, 



(7) 



with Nas the number of deposited As atoms, and 
As^' being the difference in entropy of a Ga atom 
in adsorption site i and in the Ga bulk. This dif- 
ference contains contributions from vibrational en- 
tropy as well as from configurational entropy, which 
in turn depends on the concentration of deposited Ga 
adatoms on the surface. 

A gallium atom stemming from the beam source 
typically runs through several intermediate configu- 
rations on the surface before it becomes part of the 
GaAs solid. During this process its free enthalpy gc a 
decreases gradually. Fig. g illustrates the approach of 
9Ga towards its equilibrium value in the GaAs bulk 
for the two growth scenarios discussed above. The 
quantity g Ga + As^T is shown for deposition of an 
increasing number of Ga atoms. Whenever compet- 
ing structures with or without As "capping" are ex- 
pected, we compare the free enthalpy of both. The val- 
ues shown in the figure are calculated for T = 750 K 
and p = 10 -3 Pa. Under these conditions, incorpora- 
tion of two Ga atoms in the local a structure (B3B1) 
has a lower free enthalpy than the local j3 structure 
(labelled B3B1+AS2). Similarly, larger patches of a 
structures (labelled B3B1B3BJJ have a slightly lower 
free enthalpy than the corresponding As-capped (3 
structures. Thus they will be predominantly uncov- 
ered under the considered conditions. The situation 
changes somewhat at temperatures below 700 K or 
higher As pressures, when As-covered j3 structures be- 
come equally probable. For the structures nucleated 
on the top layer, the calculations show that the AS2 
molecule binds strongly to the four-atom Ga cluster. 
Therefore the As-covered structure is preferred under 
a wide range of growth conditions where the surface 
As concentration is in equilibrium with gas phase AS2 . 
However, as can be seen from Fig. [| individual Ga 
dimers (A4 A 4 ) or Ga dimers above pairs of As dimers, 
separated by a gap (A4A 4 A 4 A 4 ), are still the most 
stable species in absolute terms. 



ft 
> 



o 



-2.9 



-3.1 



-3.3 



/ /— - n Ga (GaAs,p(As 2 )=10- Pa) 

/ / H Ga {GaAs, p(As 2 )=10~ Pa) 

/ u Ga (bulkGa) 

/ 



500 600 700 800 900 1000 
temperature [K] 

FIG. 7. Free enthalpy per Ga atom in bulk GaAs in 
thermodynamic equilibrium with As2 vapor at two pres- 
sures, compared to the free enthalpy per Ga atom in ele- 
mental Ga. 

Since the kinetics of Ga incorporation is driven 
by differences in <? Ga , a kinetic model of growth 
that includes only the Ga species should describe 

the stability of different configurations in terms of 

(i) 

the quantities g Ga . The effective parameters of such 
a model depend on temperature and arsenic back- 
ground pressure. Eq. fallows to determine them from 
first principles, apart from the yet undetermined en- 
tropic contribution AsT. In a kinetic model, all mi- 
croscopic processes are represented by their respec- 
tive rates. In addition to knowledge of the (meta- 
)stable configurations occuring during growth, a mi- 
croscopic determination of the rates also requires 
knowledge about the transition states. Within clas- 
sical rate theory, the rate for a particular transition 
is given by F = exp((<?W — g TS )/kT), with the free 
enthalpy at the transition state, g TS . It is frequently 
represented by an activation energy and a prefactor, 
r = ro exp(— Ea/HT). When rates are determined 
from calculated free enthalpies, the contributions to 
AsT from configurational entropy should be omitted, 
since they are described implicitely by the number 
of microscopic processes possible in a given situation. 
The vibrational entropies enter the rates only through 
the prefactor, but do not appear in the activation 
energy. Without a detailed calculation of vibrational 
properties, one has to rely on the assumption that 
vibrational contributions to the entropy are similar 
in different configurations. Thus they tend to cancel 
out in As and will not qualitatively change the pic- 
ture. This is equivalent to the assumption of a com- 
mon prefactor for all kinetic processes, as is frequently 
made in kinetics simulations. 
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FIG. 8. Free enthalpy per Ga atom at T = 750 K and 
p(As2) = 10 _3 Pa for two kinetic pathways marked by the 
arrows, a) for filling of the trenches (see also Fig. and 
b) for nucleation of a new layer (see also Fig. |E]) . The full 
and dashed line indicate the chemical potential of a Ga 
atom in the GaAs bulk and in elemental Ga, respectively. 



VII. CONCLUSION 

In conclusion, we have investigated the energetics of 
two kinetic pathways for homoepitaxy on GaAs(OOl) 
by means of DFT calculations. We find that two gal- 
lium adatoms on the GaAs(OOl) surface interact suffi- 
ciently strongly to form stable nuclei. One type of nu- 
clei consists of Ga atoms in adjacent three-fold coor- 
dinated sites in the trenches, while another type con- 
sists of Ga dimers forming on the top layer of arsenic. 
Since single Ga atoms in the top layer are energeti- 
cally unfavorable and constitute an excited state, the 
formation of Ga dimers in the top layer is suppressed 
at low coverages, and growth will be dominated by the 
nucleation of Ga atom pairs in the trenches. Filling of 
the trenches will proceed by further attachment of Ga 
adatoms to these nuclei. From our calculations, we ex- 
pect that islands in the new layer will preferentially 
form in those regions of the surface where locally the 
(3 reconstruction has formed, i. e. where the trenches 
have been partially filled. 

Upon further deposition of material, the islands 



start to grow. Quenched STM images have revealed 
small islands which do not yet show the (2 x 4) pattern 
[|. After passing through this metastable intermedi- 
ate state, eventually the islands will restructure and 
display the trenches characteristic for the /32 recon- 
struction. 
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